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ABSTRACT 



Planetary migration in standard models of gaseous protoplanetary disks is 
known to be very rapid (~ 10 5 years) jeopardizing the existence of planetary 
systems. We present a new mechanism for significantly slowing rapid planetary 
migration, discovered by means of radiative transfer calculations of the thermal 
structure of protoplanetary disks irradiated by their central stars. Rapid dust 
settling in a disk's dead zone - a region with very little turbulence - leaves a dusty 
wall at its outer edge. We show that the back-heating of the dead zone by this 
irradiated wall produces a positive gradient of the disk temperature which acts 
as a thermal barrier to planetary migration which persists for the disk lifetime. 
Although we analyze in detail the migration of a Super-Earth in a low mass disk 
around an M star, our findings can apply to wide variety of young planetary 
systems. We compare our findings with other potentially important stopping 
mechanisms and show that there are large parameter spaces for which dead zones 
are likely to play the most important role for reproducing the observed mass- 
period relation in longer planetary periods. 

Subject headings: accretion, accretion disks — radiative transfer — turbulence - 
planets and satellites: formation — protoplanetary disks 
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Introduction 



Extrasolar planets (ES 



s) ha ve an unexpected distribution of orbital radii around 



their host stars ( Udry et al. 

,0 



2009) - ranging from about 0.02 to 70 astronomical units 



at ion wherein 



(AU)o In particular, ESPs are observed to obey a mass-period (M-P) re 
lower mass planets end up in short period orbits around their host stars ( jUdry fc Santos 
20071 ). The predominance of very short period planets is generally thought to 
consequence of planetary migration. As an example, the tidal interaction of a planet with its 
surrounding gaseous disk excites density waves in the disk at so-called Lindblad resonances. 
These waves exert a tor que back on the planet whic h results in a net angular momentum 



transfer between them (IGoldreich fc Tremaine 



19801). P 



anets may also exchange angular 



199ll ). For locally isothermal 



momentum with the gas inside of their horseshoe region (jWardl . 
protoplanetary disks with smoothly declining distributions of disk column density and 
temperat ure with radius, the net torque generally leaves a planet spiraling inwards through 



the disk ( iTanaka et al 



20021), i.e. the 



stronger than that of the inner wake (IWard 



torque exer ted by the outer wake is marginally 



19971 ). Many calculations and simulations 



show that the migration timescale of planets in such "standard" disk models is very short 



roughl y two order s of magnitude sma 



(Myr)) (IWard 



1997 



Nelson et al. 



2000 



ler than the disk li f etime (one to ten mi 



Tanaka et al. 



2002 



D'Angelo et al. 



lion years 



20031 ). Why are 



there any planetary systems at all? 

The key to understanding the M-P relation and the survival of planetary systems is 
in how the dynamics of planetary motion is coupled to the prope rties and structure of 



the protoplanetary disks. As an example, 



Schlaufman et al. 



( 120091 ) focused on the surface 



density transition that can be produced at the location of the ic e-line, where a 



maximum can act as an accumulation point for planetesimals (IKretke fc Lin 



ocal pressure 



2007|). If it 



1 See the website: 



http:/ /exoplanet.eu 
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is assumed that type I migration is much slower than predicted in locally isothermal disk 
models, this feature could account for planets with orbital radii 0.1 - 2 AU. Obviously, a 
physical explanation for slower migration is needed. 

In this Letter, we present a new slowing mechanism of rapid type I migration -which 
may occur for planets with masses that are too low to open up a gap in their disks (massive 
planets can tidally form a gap and undergo type II migration). We show, by means of 
Monte Carlo radiative transfer simulations, th at dead zones - the dense inner disk region 



19961 ) - support a thermal barrier 



wherein turbulence cannot be readily excited (jGammid 
to migration. One of the most important consequences is that the thermal barrier could 
account for planets at larger orbital radii. In § 2, we outline our disk model and discuss 
tidal torques. In § 3, we analyze how the presence of a thermal barrier impacts the 
migration rates of low-mass planets. In § 4, we discuss potential issues for the M-P relation 
by comparing our stopping mechanism with others. 



2. Disk model & tidal torques 



Protoplanetary disks are 



star flChiang fc Goldreich 



mown to be heated by radiation from the centra 



1997 



D'Alessio et al 



1998 



Hasegawa &: Pudritz 



2010a 



hereafter, HP10). This rad iation mainly determines the thermal structure of disks 



( IKenyon fc Hartmann 



19871 ). This is because viscous heating dominates stellar irradiation 



heating only within about 1 AU for the classical T Tauri star systems (CTTSs; 



D'Alessio et al. 



19981 ) and only within about 0.1 AU for lower mass stars such as M stars 



(HP10). 



Detailed modeling of the spectral energy distributions e mitted by disks has sh own 



that s — — 1 for disks where the disk surface density £ oc r s (ID'Alessio et al 



19981 ) . It is 
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well established that the sign of a net torque exerted on planet s depends o n two central 



Ward 



19971 ). The thermal 



properties of disks, their surface density and temperature (e.g., 
structure therefore plays a critical role in controlling the direction of planetary migration. 
In disks without internal structure, the disk temperature T oc r* at the mid-plane steadily 
decreases (t < 0) and planetary migration is steadily inwards. 

The point is that disks are not simple power-law structures. The strength of 
turbulence within t hem varies con siderably, with very low levels occurring in dense regions 



called dead zones (IGammie 



(IMatsumura fc Pudritz 



onto the central star (IMatsumura et al. 



19961 ) that initially extend over roughly 10 AU in disks 



2006) and then gradu ally shrink in size as disk material is accreted 



2009 



hereafter, MPT09). Turbulence in disks is 



most likely excited by the magnetorotational instability (MRI) (iBalbus fc Hawley 



199lh 



The MRI requires good coupling between ions and magnetic fields, which is largely absent 
in the dead zone - that inner, high density region in which the ionization due to the X-rays 
from the central star and cosmic rays is suppressed. 

Dust is the dominant abso rber of stellar radiation in disks although its total mass is 100 



times smaller than that of gas (IDullemond et al. 



20091 ) . Many observations imply that it has 



a density distribution that is diff erent from the gas distribut ion, which is derived assuming 



vertical hydrostatic eq uilibrium (IKenvon fc Hartmann 



of its size distribution (IDullemond &: Dominik 



19871 ). Dust settling, a consequence 



2004at ). is ubiquitous in protoplanetary disks 



around any young star (HP10, references herein). The dust scale height depends upon the 
amplitude of disk turbulence which keeps it suspended in the g ravitational field of the disk 



(which is determined by the central star) (iDubrulle et al.lll995l ). 



We demonstrated in HP10 that dust settling in dead zones results in the appearance 
of a limited region in which the disk temperature can actually increase with radius - a 
radial temperature inversion. In this Letter, we adopt the disk model developed in HP 10 
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(see Table. [T]) and focus our computations on M dwarf systems, such as the recently 
discovered Super-Ear th (~ 5M e ; M e is the Earth's mass) with an orbital radius of 2 AU 



( IBeaulieu et al 



20061 ). We refer the readers to HP10 for the detail. The low mass of disks 
around M stars allows much more comprehensive Monte Carlo simulations to be performed, 
but our analysis in principle applies to any protoplanetary disk. 



We adopt the torque formula (IWard 



1997 



Menou fe Goodman 



2004 



Jang-Condell fc Sasselov 



20051 ) in which only the Lindblad torque is considered. This is because the corotation torque 
is readily saturated (i.e. is canceled out) in dead zones in our radiatively heated disk model. 
We compare the libration timescale (ie the timescale for gas to complete an orbit in the horse- 
shoe region), tuj, ~ 8nr p /3Q p x s , with the viscous timescale, jj/3z/, where the half- 



widt 



i of the horseshoe region is x s /r p = l.Q8( M„r p / 'M^hp) 1 ^ 2 (Paardekooper fc Papaloizou 



20091 ). the kinematic viscosity is v = ah 2 flKe P ( jShakura fc Sunyaev 



19731 ). and the disk 



scale height is h. Our numerical results give h p /r p ~ 0.05, so that the critical value of 
turbulence, a cr u, below which the corotation torque is saturated (and therefore negligible), 
is 

/ U V' 2 / U \ - ' 2 /h p /r p y 7/2 



a 



crit 



0.01 ' p 



5M ff 



O.IMp 



0.05 



Since the dead zone has a = 10 -5 , we can safely neglect the corotation torque for 
M p > 0.5M e in the dead zone of our disk model. In addition, we confirmed that, in the 
active region where the corotation torque is generally unsatura ted, both Lindblad a nd 



2009|). 



corotation torques result in inward migration in our disk model (jPaardekooper et al. 
Thus, exclusion of the corotation torque in the active region does not affect our findings 
in our disk model. We note that corotation torque may be unsaturated in dead zones for 
sufficiently small planetary masses, but the exact limit will depends on knowing the disk 
scale height that is established by the host star. Furthermore, our torque formula takes into 
accou nt the effects of vertical disk th ickness by diluting the gravitational force of a planet 



by z (jJang-Condell fc Sasselov 



20051 ) . 
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Results 



We performed numerical simulations of the thermal str ucture of 2D disks by solving 



2004b 



the radiative transfer equation with a Monte Carlo method (IDullemond &: Dominik 
HP10). We included the effects of vertical hydrostatic balance, dust settling, a dead zone, 
and the gravitational field of a planet embedded in the disk. The tidal torque is calculated 
as described in § 2, and incorporates our numerical data, in order to calculate the migration 
time. 

Figure [1] presents the thermal and density structure of a disk with a 5M e planet placed 
at 8 AU. The top and bottom panels show the dust and gas densities by color, respectively. 
Since we define the disk temperature as the mass-averaged temperature of dust, both panels 
show the identical temperature structure which is represented by contours. The thick line 
denotes the Hill radius r# = r p (M p /3M^) 1 ^ 3 . In this Letter, we adopt, without loss of 
generality, a dead zone which is 6 AU in size. Dead zones enhance dust settling because of 
the low turbulent amplitude there. Since disks have inner dead, and outer active regions for 
turbulence, the transition of the density distribution of dust occurs at the outer edge of the 
dead zone. This leaves a marked step in the dust scale height behind - in effect a wall of 
dust. The additional stellar e nergy absorbed at the wall is distributed by radiative diffusion 



( jHasegawa fc Pudritz 



2010bl ) since the optical depth at this region is high. The resulting 
radial " thermal inversion" - ie a region of increasing temperature with increasing disk radius 
- is shown in Figure HJ An analytic fit to our data shows that this back-heated region in 
the dead zone has a positive temperature gradient described by a power-law T oc r*' with 
f > 3/2. 

We show in Figure [3] (top) that this radial thermal inversion causes the migration rate 
to be positive (the migration time becomes negative - bottom panel), - i.e. planets migrate 
outward in the region with the positive temperature gradient. The physical explanation of 
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this behavior is that the increasing function of disk temperature changes the disk's pressure 
distribution which in turn causes the position of the outer Lindblad resonances to be further 



19931 ). This results in outer torques that 



from the planet than the inner ones ( Artymowica 
are much weaker than the inner oneso Figure [3] (bottom) also shows that the planets very 
slowly enter the region of torque reversal - as seen by the strong positive " spike" in the 
migration time. These regions correspond to radii at which dT/dr ~ (see Figure [2]) at the 
inner and outer resonances, making the torque difference between them very small. 

We emphasize that the positive temperature gradient arising from the wall-like dust 
structure is achieved for the case of a finite transition region, Ar < 10/t, in the value of 
turbulent a (HP10). Although, for simplicity, we adopt a sharp spatial transition from the 
active to the dead zone in this Letter (Ar = 0), our above results are valid for the case of 
Ar ~ h, since the positions of Lindblad resonances are typically offset from the planets 



by 2h/3 (lArtymowiczl Il993l ) . In addition, it is interesting that the migration timescale of 
the M star system is similar to that of CTTS (~ 10 4 — 10 5 years) for the other two cases 
(well mixed, dust settling). This is becaus e the tidal torque is s c aled by £(/i/r)~ 2 . A more 



detailed discussion of them is presented in 



Hasegawa fc Pudritzl (l2010bl ). 



4. Discussion 

4.1. A thermal vs a density barrier at outer dead zone radius 

We find that the dusty wall produces a radial temperature inversion that is a thermal 
barrier for rapid type I planetary migration. Whereas the torque balance in the well 
coupled active zone forces planets to migrate inward, once they encounter the radial thermal 



2 We will sho w in a forthcoming paper th at the Lindblad torque becomes negative when 
s-t/2< -7/4 flHaseeawa fc Pudritzl boiObh . 
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inversion region, the torque balance reverses, and they move out of the region. Thus, 
planets are trapped there if they originally migrate from the active region beyond the dead 
zones, or even if they formed close to the outer edge of the dead zone. 

The astrophysical implications of this result are very important since we have shown 
that dusty protoplanetary disks with dead zones possess an innate mechanism for strongly 
slowing planetary migration within them, provided that the corotation torque is saturated. 
While such a thermal barrier exists for any size of dead zone (HP10), its effectiveness is 
probably most important for lower mass disks, as we now show. 

The density structure of disks evolves with time due to viscous evolution. MPT09 
found that the difference of a between the active and dead regions produces a steep density 
gradient at the boundary and the location of their jump moves inward with time over the 
long (~ 10 Myr) viscous timescale of the dead zone. This density gradient at the outer dead 
zone boundary also plays an important role in slowing down or stopping type I migration, 
provided that planets migrate from larger disk radii. The inner torques become larger than 
the outer ones in the density gradient region, resulting in the reflection of migrating planets 
off the density gradient (MPT09). 

The relative importance of these two dead zone mechanisms is controlled by the ratio 
of dust settling r set w E/ \/2ii p d aVl Kep and the viscous r vis timescales, where p d is the 
bulk density of dust and a is the grain size of dust. We find that the critical condition 
Tset/Tvis > 1 for dominance of the density vs thermal barriers presented by a dead zone is 

S (?) " 25 (lF*) (li^) (lrn^i) = ° cnt - (2) 
This implies that a density barrier is dominant for sufficiently large values of £ and h/r. 
For disks around CTT Ss with a typical dead zon e size (~ 10 AU), £(/i/r) 2 ~ 300 g 



cm" 2 x (0.4) 2 w 2C crit (jChiang & Goldreich 



19971 ). which indicates that a density barrier is 



dominant. For disks around M stars with a typical dead zone size (~ 5 AU), £(/i/r) 2 ~ 20 
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g cm" 2 x (0.05) 2 = 2 x lQ- 3 C crit flScholz et al. 



20071 ). which implies that a thermal barrier 



is dominant. Generally, we find that a thermal barrier becomes weaker in the late stages 
because of accretion which reduces the density at the outer edge of the dead zones. 



4.2. Comparisons with other possible stopping mechanisms and potential 

effects on the M-P relation 



It is wel 



central star ( iLin et al. 



known that ti dal interaction and angular momentum exchange with the 
19961 ) and th e creation of a hole in the inner part of disks by the 



presence of a stellar magnetosphere (IShu et al. 



1994 ; 



Lin et al. 



19961 ) cannot reproduce the 



whole of the observed M-P relation. This is because such barriers become important only 



for planets approaching very close to the cent ral star. Most 



orbital radii from about 0.02 AU to 10 AU ( jUdry fc Santos 



SPs, however, have observed 



20071 ). Thus, while these 



barriers may play a role in piling up ESPs in the vicinity of the star, they have difficulty in 
predicting planets with larger orbital radii. 

Stochastic migration provides another mechanism for contr olling planetary migration. 



It arises when disks undergo magnetohydrodynamic tur bulence (Nelson fc 



Laughlin et al. 



20041 ). which is the outcome of the MRI ( IBalbus & Hawlev 



Pa paloizou 



2004 



1991) . Stochast ic 



2006). 



torques tend to reduce the timescale for planetary survival, however ( 1 Johnson et al. 
In a few exceptional cases, planets can diffuse out to large disk radii where they can survive 
longer. Planets within the dead zones cannot perform random walks because turbulent 
torques are reduced by about two orders of magnitude th ere and consequently would 



undergo steady inward type I migration ( lOishi et al. 



20071 ). Thus, stochastic migration is 



unlikely to be the main barrier to rapid type I migration. 



Is our thermal barrier sufficiently wide to stop planets scattered into the dead zone by 
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stochastic effects? We consider a characteristic length scale for turbulent diffusion defined 
by Ar turb = y/vT c = h^/aT r /VL~ l , where r r is th e correlation time of turbulence. We adopt 



a value r c = 0.5f2 1 ( INelson fc Papaloizou 



20041 ) and find that Ar turb » 0.02 AU at r = 6 



AU. This is shorter than the width of the thermal barrier (~ 2 AU). Therefore, stochastic 
motions due to turbulence are unlikely to affect the migration stopping mechanism at the 
thermal barrier. 

Corotation torques are also potentially important for slowing or stopping planets. In 
both barotropic and adiabatic disks, corotation torque associated with a radial vortensity 
grad ient may work as a barrier around the re gion of inner edge of the dead zone (~ 0.01-0.1 



AU; 



Fromang et al. 



2002 



Masset et al. 



2006 



, MPT09) because it becomes po sitive due to 



a pos itive surface density gradient, resulting in outward planetary migration (iMasset et al 



20061 ) . However, the location of the barrier is almost constant with time because stellar 
irradiation controls it, and consequently this barrier is important only for planets in the 
vicinity of the star. In adiabatic disks, corotation torque associated with a radial entropy 
gradient may also act as a barrier because it becomes la rge and positive around the regi on 



with a large (in magnitude), negative entropy gradient ( 



Baruteau fc Masset 



2008 



Paardekooper fc Mellema 



2008 



aardekooper fc Mellema 



2006 



Paardekooper fc Papaloizou 



20081 ). Thus, corotation torque may be important to the M-P relation in certain regimes of 
planetary mass and disk heating, as noted above. 

Our mechanism has a movable barrier that shrinks from larger disk radii on the 10 
Myr (for CTTSs) viscous timescale of the dead zone. As a concrete example, this shrinkage 
of the dead zone cou ld explain the recent detected Super-Earth at 2 AU around an M star 



(Be aulieu et al. 



20061 ). This is because a Super-Earth is likely to be the most massive planet 
that would surely form in a low mass disk and is therefore most likely to have been formed 
beyond the outer dead zone radius and left behind as the dead zone shrinks away. 
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We conclude that the stopping mechanisms arising from dead zones - via thermal and 
density gradients - are important barriers to rapid planetary migration. We have shown 
that the thermal barrier that arises from disk irradiation by a heated dust wall is robust 
and may be most important in the earlier phases of disk evolution and for the evolution 
of low mass systems as found around M stars as an example. Dead zones may provide 
the most significant slowing mechanism of type I migration that is needed to explain the 
longer period of population of the M-P relation, which will be checked in future population 
synthesis models. 
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by McMaster University, as well as by Graduate Fellowships from SHARCNET and the 
Canadian Astrobiology Training Program (CATP). REP is supported by Discovery Grant 
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Table 1: Summary of parameters fc symbols 



Symbol 


Meaning 


Value 


M* 


Stellar mass 


0.1 M 


R* 


Stellar radius 


O.4R 




Stellar effective temperature 


2850 K 


£ 


The surface density of (gas + dust) 


oc r _1 


M d 


The total disc mass (gas + dust) 


4.5 x 10~ 3 M 




gas-to-dust ratio 


100 


a 


Parameter for turbulence (active/dead) 


10- 2 /10- 5 


M p 


Planetary mass 


5 M e 



M is one solar mass, and R Q is one solar radius. 
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Fig. 1. — The density and temperature structures of disks with a 5M ffi mass planet located 
at 8 AU. For both panels, the density is denoted by colors [g cm -3 ] as shown in colorbar, and 
the disk temperature is denoted by contours [K] (both panels show the identical temperature 
since the disk temperature is defined by taking the mass-average of the dust temperatures). 
The size of a dead zone is 6 AU. Top: the dust density. Bottom: the gas density. The thick 
black circle denotes the Hill sphere r# = r p (Mp/3M») 1,/3 , where r p and M v is the position 
and mass of a planet, respectively, and M* is the stellar mass. A wall-like structure appears 
at the boundary between the active and dead zones in the dust distribution while it is not 
in the gas. 



-18 - 



1000 



10 

















r 3/7 (Chiong & Go 


dreich model) 








W/ Dust settling & 


a dead zone 








r V2 (Analytical fit) 
































































< 











2 4 6 8 10 

r (AU) 



Fig. 2. — The disk temperature in the disk mid-plane. The solid line is for the case of dea d 



1993), 



zone, the dashed line is the analytical model of disk temperature ( iChiang fc Goldreichl l 
and the dashed-dot line is for analytical fit to a positive temperature gradient. The size of 
a dead zone, which is 6 AU, is indicated by the vertical solid line. 
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Fig. 3. — Migration rate and time for a 5 mass M e planet at various orbital radii on the top 
and bottom panels, respectively. For both panels, the solid line denotes the case of a dead 
zone which is 6 AU in size whose position is indicated by the vertical solid line, the red dashed 
line is for the well-mixed case, and the blue dashed-dot line is for the dust settling case. The 
negative migration time region in the bottom panel (outward migration) corresponds to the 
region with a positive temperature gradient (see Figure [2]). 



